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[i] We evaluate the sensitivity of tropospheric OH. 0 3 , and 0 3 precursors to 
photochemical effects of aerosols not usually included in global models: ( 1 ) aerosol 
scattering and absorption of ultraviolet radiation and (2) reactive uptake of H0 2 , N0 2 , and 
NOj. Our approach is to couple a global 3-D model of tropospheric chemistry (GEOS- 
CHEM) with aerosol fields from a global 3-D aerosol model (GOCART). Reactive uptake 
by aerosols is computed using reaction probabilities from a recent review (y H02 = 0.2, 

*Yno 2 = 1CT 4 . Yxo 3 = 10-'). Aerosols decrease the 0 3 — O('D) photolysis frequency by 
5-20% at the surface throughout the Northern Hemisphere (largely due to mineral dust) 
and by a factor of 2 in biomass burning regions (largely due to black carbon). Aerosol 
uptake of H0 2 accounts for 10-40% of total HO x radical (= OH + peroxy) loss in the 
boundary' layer over polluted continental regions (largely due to sulfate and organic 
carbon) and for more than 70% over tropical biomass burning regions (largely due to 
organic carbon). Uptake ofN0 2 and N0 3 accounts for 10-20% of total HN0 3 production 
over biomass burning regions and less elsewhere. Annual mean OH concentrations 
decrease by 9% globally and by 5-35% in the boundary' layer over the Northern 
Hemisphere. Simulated CO increases by 5-15 ppbv in the remote Northern Hemisphere, 
improving agreement with observations. Simulated boundary layer 0 3 decreases by 15- 
45 ppbv over India during the biomass burning season in March and by 5-9 ppbv over 
northern Europe in August, again improving comparison with observations. We find 
that particulate matter controls would increase surface 0 3 over Europe and other industrial 
regions. INDEX TERMS: 0305 Atmospheric Composition and Structure: Aerosols and particles (0345. 

4801): 0345 Annospheric Composition and Structure: Pollution — urban and regional (0305): 0360 
Atmospheric Composition and Structure: Transmission and scattenng of radiation; 0365 Atmospheric 
Composition and Structure: Troposphere — composition and chemistry; KEYWORDS: aerosols, ozone, 
heterogeneous chemistry, radiation, reactive uptake. OH 
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1. Introduction 

[:] Understanding the radiative and heterogeneous 
chemical effects of aerosols on OH and 0 3 , the main 
tropospheric oxidants, is a major issue for models of the 
oxidizing power of the atmosphere, radiative forcing, and 
surface air quality. The importance of N 2 0 5 hydrolysis in 
aerosols is well established [He ikes and Thompson. 1983: 
Dentener and Crutzen , 1993; Tie et ah. 2001: Liao et ah. 
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2003] and is included in all current-generation global 
models of tropospheric chemistry. Model calculations by 
Dentener and Crutzen [1993] indicate that this reaction 
reduces global tropospheric 0 3 and OH by 9%. Other 
aerosol effects have received far less study and are the 
focus of this paper. 

[ 3 ] Scattering and absorption of ultraviolet (UV) radia- 
tion by aerosols modifies photolysis rates [He and Carmi- 
chael, 1999; Liao et ah. 1999], Reactive heterogeneous 
uptake of H0 2 , N0 2 , and N0 3 by aerosols may be impor- 
tant as well [Jacob, 2000], Including these photochemical 
effects of aerosols in global models could help resolve some 
long-standing discrepancies with measurements, such as the 
underestimate of the CO background in the Northern Hemi- 
sphere [Hauglustaine et ah, 1998; Lawrence et ah, 1999; 
Bergamaschi et ah, 2000. 2001: Bey et ah, 2001a] and the 
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overestimate of 0 3 over India [Lai and Lawrence, 2001; 
Martin et al., 2002b]. 

[4] Recent development of global simulations and obser- 
vational databases for the major aerosol types (sulfate, black 
carbon, organic carbon, sea salt, and mineral dust) [Tegen et 
al., 1997; Takemura et al., 2000; Glum et al., 2001; Chin et 
al., 2002; Chuang et al., 2002; Liao et al., 2003] facilitates 
a more detailed investigation of the photochemical effects of 
aerosols on tropospheric 0 3 and OH. We present here such 
an investigation using a global 3-D model of tropospheric 
chemistry (GEOS-CHEM). Section 2 describes the imple- 
mentation of aerosol photochemistry in GEOS-CHEM and 
section 3 examines its implications. Section 4 compares the 
results with observations for species and regions where the 
effects are particularly large. Global budgets of 0 3 and OH 
are analyzed in section 5. 

2. GEOS-CHEM Model 

2.1. General Formulation 

[ 5 ] The GEOS-CHEM model was initially described by 
Bey et al. [2001a] and subsequent improvements are 
described by Martin et al. [2002b]. We use here GEOS- 
CHEM version 4.26 (http://www-as.harvard.edu/chemistry/ 
trop/geos/) with improved treatment of aerosol photochem- 
ical effects and a few other updates as described below. The 
model is driven by assimilated meteorological data from the 
Goddard Earth Observing System (GEOS) of the NASA 
Data Assimilation Office (DAO) [Schubert et al., 1993]. 
The meteorological data include 3-D fields updated every 3 
hours for surface fluxes and mixing depths and every 6 
hours for other relevant variables. We use for this study the 
GEOS data for 1996-1997, available with a resolution of 2° 
latitude by 2.5° longitude and 46 sigma levels in the vertical 
extending up to 0.1 hPa. For computational expedience we 
degrade the horizontal resolution to 4° latitude by 5° 
longitude and merge the vertical levels above the lower 
stratosphere, retaining a total of 26. There are five levels 
located below 2 km altitude (for a column based at sea level, 
these levels are centered at 50 m, 250 m, 600 m, 1100 m, 
and 1700 m). 

[ 5 ] The GEOS-CHEM model includes a detailed descrip- 
tion of tropospheric 0 3 -NO x -hydrocarbon chemistry. It 
solves the chemical evolution of about 120 species with a 
Gear solver [. Jacobson and Titrco, 1994] and transports 24 
tracers. Photolysis frequencies are computed using the Fast- 
J radiative transfer algorithm [ Wild et al., 2000] which 
includes Rayleigh scattering as well as Mie scattering by 
clouds and aerosols. Extensive evaluations of the GEOS- 
CHEM chemical fields with observations of tropospheric 
0 3 and its precursors in different regions of the world have 
been presented in a number of papers [e.g., Bey et al., 
2001a, 2001b; Li et al., 2001, 2002a, 2002b; Palmer et al., 
2001, 2003; Chandra et al., 2002; Fiore et al., 2002; Liu et 
al., 2003; Martin et al., 2002a, 2002b]. Observed monthly 
mean 0 3 concentrations in different regions of the tropo- 
sphere are usually reproduced to within 10 ppbv with no 
global bias, but there are some regional problems, most 
prominently a 20-30 ppbv positive bias over India [Martin 
et al., 2002b]. The seasonal and interannual variability of 
CO is well captured, but there is a general negative bias 
of 10-20 ppbv in the Northern Hemisphere. Consideration 


of aerosol photochemical effects alleviates these discrep- 
ancies, as discussed in section 4. 

[7] The model version used here includes several updates 
to the version presented by Martin et al. [2002b], in 
addition to the aerosol effects described below. The most 
important are updated reaction rates for organic peroxy 
radicals [Tyndall et al., 2001], updated CO emissions 
(B.N. Duncan et al., Model study of the variability and 
trends of carbon monoxide (1988-1997), 1, Model formu- 
lation, evaluation, and sensitivity, manuscript in preparation, 
2002), and updated acetone emissions [Jacob et al., 2002], 
Also, surface emissions of trace gases are now distributed 
vertically in the unstable mixed layer (which can vary from 
less than 100 m at night to over 2000 m during the day) to 
minimize numerical artifacts from operator splitting (1-hour 
time steps between chemistry and transport). 

2.2. Aerosol Radiative and Chemical Effects 

[x] Previous versions of GEOS-CHEM have included 
heterogeneous reactions of HO 2 , N0 2 , N0 3 , and N 2 Oj in 
sulfate aerosols [Bey et al., 2001a], and more recently these 
same reactions as well as radiative effects have been 
included for mineral dust aerosols [Martin et al., 2002b]. 
We include here a more comprehensive treatment of the 
photochemical effects of aerosols. To this end we use global 
3-D monthly mean fields for 1996-1997 of dry aerosol 
mass from the GOCART model [Chin et al., 2000a, 2000b. 
2002; Ginoux et al., 2001] which uses the same GEOS 
meteorological fields and transport algorithms as GEOS- 
CHEM. The GOCART data include sulfate, size-resolved 
mineral dust and sea salt, hydrophobic and hydrophilic 
black carbon, as well as hydrophobic and hydrophilic 
organic carbon. Nitrate aerosol is not included. 

[4] As described in detail elsewhere [Chin et al., 2000a, 
2002; Ginoux et al., 2001], the GOCART model sulfate 
simulation includes oxidation of S0 2 emitted from indus- 
trial activity, biomass burning, biofuel, volcanic eruptions, 
and formed from dimethylsulfide (DMS) released from the 
ocean, as well as a small direct sulfate source. Oxidation 
rates in the GOCART model are calculated with offline 
fields of OH and aqueous H 2 0 2 . Organic and black carbon 
are emitted directly from biomass burning and industrial 
activity, with 80% of black carbon and 50% of organic 
carbon emitted as hydrophobic, the rest being hydrophilic. 
Organic carbon is also fonned from volatile organic com- 
pounds. Hydrophobic carbonaceous aerosols age to become 
hydrophilic with an e-folding time of 1 .2 days. Mineral dust 
emission is calculated as a function of the local probability 
of bare sediments, surface wind speed, and surface wetness. 
Sea-salt emission is calculated as a function of surface wind 
speed. All aerosols are subject to dry deposition and hydro- 
philic aerosols to wet deposition. 

[ 10 ] Simulation of the radiative and chemical effects of 
aerosols requires information on the relative humidity (RH) 
dependent aerosol size distribution. We obtain this informa- 
tion from Ginoux et al. [200 1 ] for mineral dust, and from 
the Global Aerosol Data Set (GADS) [Kopke et al., 1997] 
and Chin et al. [2002] for other aerosols, and apply it to 
local RH from the GEOS fields. Table 1 shows the hygro- 
scopic growth factors and Table 2 shows the effective radius 
r e at 70% humidity. The effective radius is defined as 
j,.'Kr i f(r)dr/f r Tu- 2 f(r)dr, where /(r) is the fraction of par- 


MARTIN ET AL.: PHOTOCHEMICAL EFFECTS OF AEROSOLS 


ACH 6 - 3 


Table 1. Hygroscopic Growth Factors for Aerosols 3 





RH. % 



0 

50 

70 

80 

90 

Sulfate 

1.0 

1.4 

1.5 

1.6 

1.9 

Black carbon h 

1.0 

1.0 

1.0 

1.2 

1.4 

Organic carbon h 

1.0 

1.2 

1.4 

1.5 

1.7 

Sea salt 

1.0 

1.6 

1.8 

2.0 

2.4 

Mineral dust c 

1.0 

1.0 

1.0 

1.0 

1.0 


“The hygroscopic growth factor is the ratio of the wet to the dry effective 
radius. 

b These are values for hydrophilic aerosols. There is no hygroscopic 
growth for hydrophobic carbon aerosols. 

c Dust particles are sufficiently large, and their hygroscopicity sufficiently 
low. that hygroscopic growth is considered negligible [Li-Jones et al. 
1998], 


tides w ith radius between r and r + dr [Hansen and Travis , 
1974]. We cap aerosol growth above 90% RH: hygroscopic 
growth is very nonlinear at higher RH and cannot be 
adequately resolved in the GEOS fields because of sub- 
grid-scale v ariability, especially associated with clouds. 

[i l ] Calculation of the radiative effects of aerosols also 
requires information on the w'avelength-dependent complex 
refractive index. We obtain this information from Patterson 
et al. [1977] for mineral dust, and for other aerosols as a 
function of RH from GADS. We assume an external 
mixture of the different aerosol types. Posfai et al. [1999] 
found that soot and sulfate aerosols arc often internally 
mixed: Liao et al. [1999] calculated that the absorbing effect 
of aerosols on photolysis frequencies is 5-15% larger for 
internally versus externally mixed sulfate-soot aerosols. A 
more physically realistic treatment of mixed sulfate-soot 
aerosols as a black carbon core surrounded by sulfate would 
have an intermediate effect [Jacobson. 2000]. 

[ 12 ] We use a Mie algorithm [de Rooij and van der Stap, 
1984; Mishchenko et al., 1999] to calculate for each aerosol 
type the single scattering albedo, the extinction efficiency, 
and the first eight terms in the Legendre expansion of the 
phase function, for input to the Fast-J radiative transfer code 
of Wild et al. [2000]. We tabulate the optical properties at 
the RH values in Table 1. and interpolate in the simulation. 
Table 2 shows the resulting optical properties at 0.4 pm and 
70% RH. Sulfate and sea salt are purely scattering, black 
carbon is strongly absorbing, and organic carbon is moder- 
ately absorbing at UV and visible wavelengths. Organic 
carbon aerosols could be highly absorbing in the UV if they 
contain chromophore functional groups [Jacobson. 1999], 
but this is not considered here. Mineral dust is strongly 
absorbing in the UV and moderately absorbing at visible 
wav elengths. The extinction efficiencies of black carbon, 
organic carbon, and sulfate are about 20-60% larger at 300 
nrn than at 400 nm. The extinction efficiencies of super- 
micron mineral dust and sea-salt particles exhibit little 
wavelength sensitivity at UV and visible wavelengths. 

[ 13 ] We calculate the aerosol optical depth locally in the 
model from the mass concentration, extinction efficiency, 
and particle mass density for each particle type and wave- 
length. Figure 1 shows the modeled aerosol optical depth at 
400 nm during March and August. Extensive evaluation of 
the GOCART fields w'ith in situ, ground based, and satellite 
observations have been presented in a number of papers 
[Chin et al., 2000b. 2002: Ginoux et al., 2001 : Penner et al.. 


2001; S. Kinne et al.. Monthly averages of aerosol proper- 
ties: A global comparison among models, satellite data and 
AERONET ground data, submitted to Journal of Geophys- 
ical Research, 2002]. The GOCART aerosol fields are 
generally consistent w'ith observations to within a factor 
of 2. Mineral dust and sulfate dominate the background 
optical depth in the Northern Hemisphere. Mineral dust 
optical depth reaches a seasonal maximum during boreal 
summer, associated with surface cyclones that erode par- 
ticles from topographic depressions in desert regions and 
with intense solar heating that transports particles into the 
free troposphere through dry convection [Ginoux et al.. 
2001]. The mineral dust optical depth of greater than 0.5 
over much of northern Africa during summer is consistent 
with measurements from AERONET [Holben et al., 2001] 
and retrievals from TOMS [Torres et al.. 2002], as shown 
by Ginoux et al. [2001] and Chin et al. [2002], Mineral dust 
over the tropical Pacific is biased high by a factor of 2-5 
[Ginoux et al.. 2001], Sulfate optical depth is maximum 
during summer due to higher S0 2 oxidation rates, and is 
most important over the eastern United States, eastern Asia, 
and Europe [Chin et ah. 2000a. 2000b]. Modeled sulfate 
concentrations over the eastern United States are a factor of 
1.5 lower than observations from the Eulerian Model 
Evaluation Field Study (EMEFS) [Chin et al., 2000b]. 
Black and organic carbon aerosols from biomass burning 
contribute to seasonal enhancements in optical depth over 
South Asia during March and southern Africa during 
August [Kuhlbusch et al., 1996: Chowdhury et al.. 2001; 
Holben et al, 2001; Leon et al., 2001]. The GOCART 
model optical depths are generally within 0.3 of AERONET 
measurements over African biomass burning regions, but 
are about 2-3 times lower than AERONET measurements 
over the biomass burning region of South America [Chin et 
al.. 2002]. Black carbon is also notable over northern 
Europe (diesel engines) and East Asia (biofuels). Sea-salt 


Table 2. Aerosol Optical Properties at 0.4 pm Wavelength (RH = 

70%) a 


Aerosol Type 

Density', g cm 

r e . pm 

Refractive Index 

-'-O 

0 

Sulfate 

1.7 

0.24 

1.36-10% 

1.00 

1.95 

Black carbon 

1.0 

0.04 

1.75-0.46/ 

0.27 

0.76 

Organic carbon 

1.8 

0.10 

1.42-0.002/ 

0.98 

0.74 

Sea salt b 

Accumulation 

■> ”> 

1.3 

1.37-10% 

1.00 

2.43 

Coarse 

2 2 

10.1 

1.37-10% 

1.00 

2.09 

Mineral dust b 

0.1-0.18 pm 

2.5 

0.15 

1.58— 0.014i 

0.94 

2.48 

0.18-0.3 pm 

2.5 

0.25 

1 .58—0.0 1 4i 

0.91 

3.10 

0.3-0. 6 pm 

2.5 

0.4 

1.58-0.0141 

0.86 

2.82 

0.6-1 pm 

2.5 

0.8 

1 .58 — 0.014i 

0.76 

"■t 
r i 

1-1.8 pm 

2.65 

1.5 

1.58— 0.014i 

0.68 

2 *>7 

1.8-3 pm 

2.65 

2.5 

1 .58 0.0 1 4i 

0.62 

2.19 

3-6 pm 

2.65 

4.0 

1 .58— 0.014i 

0.58 

2.14 


“Effective radius r e . single scattering albedo and extinction efficiency 
Q are calculated using refractive index and lognormal size distribution data 
available from the Global Aerosol Data Set (GADS) [Kopke et at.. 1997] 
and Chin et al. [2002] with the following exceptions. We use a geometric 
standard deviation of 2.0 for all aerosol types except mineral dust. The 
calculations for mineral dust use complex refractive indices from Patterson 
et al [1977]. and assume a standard gamma particle size distribution with 
effective variance of 0.2. 

b Size ranges were resolved in the GOCART simulation and used as input 
to GEOS-CHEM. 
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Figure 1 . Modeled monthly mean optical depth at 400 nm of different aerosol types for March 1997 
(left column) and August 1997 (right column). Values are calculated from mass concentration fields from 
Ginoux et al. [2001] for mineral dust and from Chin et al. [2002] for the other aerosol types. See color 
version of this figure at back of this issue. 


optical depths are highest near ±60° where strong surface 
winds increase emissions [ Gong et al 1997]. 

[ 14 ] We implement heterogeneous chemistry on aerosol 
surfaces in GEOS-CHEM following a standard reaction 
probability formulation. We neglect heterogeneous chemis- 
try in clouds which would be difficult to constrain in our 
coarse resolution model; its effect on O 3 and OH is limited 
in any case by the small fraction of atmospheric volume 
occupied by clouds [Liang and Jacob, 1997], As pointed 
out by Jacob [2000] it is the atmospheric volume fraction 
occupied by clouds, rather than the frequency of cloud 
processing, that limits the effect of cloud chemistry on the 
global budgets of 0 3 and OH. 


[15] In the reaction probability formulation, the rate 
constant k for chemical loss of a gas with mean molecular 
speed v and gas-phase molecular diffusion coefficient D g on 
an aerosol of particle radius a is given by 



where y is the reaction probability, i.e., the probability that a 
molecule impacting the aerosol surface undergoes reaction 
[ Ravishankara , 1997], and A is the aerosol surface area 
per unit volume of air. We calculate D g as a function of 
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Figure 2. Modeled monthly mean total aerosol surface area 
in the lower troposphere (~600 m above the surface) for 
March and August 1997. Values are calculated from mass 
concentration fields from Ginoux et al. [200 1 ] for mineral 
dust and from Chin et al. [2002] for the other aerosol types. 
See color version of this figure at back of this issue. 


molecular weight, temperature, and air density following 
Dentener [1993]. We include the four heterogeneous 
reactions H0 2 — * 0.5 H 2 0 2 . N0 2 — > 0.5 HONO + 0.5 
HN0 3 , NO : , — HNO;. and N 2 0 5 — 2 HNO; with reaction 
probabilities of 0.2. 10~ 4 , 10~ 3 , and 0.1. respectively, as 
recommended by Jacob [2000]. The products have long 
lifetimes relative to the timescales for gas-particle exchange 
and are therefore released into the gas phase. Uptake of N0 2 
on dry dust surfaces may be less (10 -6 -10 -4 ) [Grassian. 

2001] , but the phase of mineral dust surfaces in the 
atmosphere is unclear [Martin, 2000], For the range of r e 
and y values used here, the chemical rate constant k may be 
limited either by free molecular collision (4 \~y term in 
equation (1)) or by diffusion (a/D g term). Figure 2 shows 
the aerosol surface area in the lower troposphere (~600 m). 
The spatial distribution is similar to that of optical depth but 
features stronger maxima in regions where submicron 
aerosols dominate. 

[i6] We do not include reactive uptake of 0 3 by aerosols, 
for which Jacob [2000] cites no conclusive laboratory or 
field evidence. Dentener et al. [1996] postulated direct 
uptake of O; by dust with an assumed y 0 ; of 10 *—10 ' 
based on analogy with measured 0 : , deposition to bare soil. 
However. O; deposition to soil most likely involves reaction 
with unsaturated organics. Laboratory measurements show 
that Oj uptake by organic aerosol surfaces is rapidly 
quenched as surface reaction sites are oxidized [Moise 
and Rudich, 2000]. Recent laboratory measurements on 
dry mineral dust surfaces reveal 703 ~ 10 - ' with no 
apparent surface saturation, but with a factor of 2-3 
decrease over several hours of exposure [Michel et al., 

2002] . A model study of an African dust layer over the 
Atlantic by De Reus et al. J2000] indicates that reactive 
uptake of O; (7 = 5 x 10 ^1 would reduce simulated 0 3 
concentrations in the layer by 1 1 ppbv, improving agree- 


ment with observations. Further work is evidently needed 
on this issue. 

[r] We also do not include reactive uptake of CH 3 0 2 . 
HCHO. HNO3, or H 2 0 2 by aerosols. Current evidence is 
that uptake of CH 3 0 2 and HCHO is too slow to be of 
atmospheric significance [Jacob, 2000], Uptake of HNO; 
by alkaline dust particles is important for the partitioning of 
HNO; between the aerosol and gas phase, but should have 
little impact on oxidant concentrations since both HNO;(g) 
and nitrate aerosols are removed by wet deposition. De Reus 
et al. [2000] calculated that uptake of HNO; had a large role 
only if the HNO; in the aerosol phase was not permitted to 
regenerate NO x through photolysis. Uptake of H 2 0 2 , if 
determined by reaction with S0 2 . is significant only in 
clouds and has little overall effect on HO x [Tie et al., 2001 ]. 

[is] We use a March 1996 to August 1997 simulation 
period. The first 6 months are used to achieve proper 
initialization and we present results for September 1996 to 
August 1997. Biomass burning emissions for that period are 
specified from satellite observations [Duncan et al., 2003] 
(the large fires in Oceania associated with El Nino have not 
yet begun). The standard simulation is as described above. 
A sensitivity simulation is conducted that excludes both the 
radiative effects of aerosols and the heterogeneous uptake of 
H0 2 . N0 2 , and NO; (but retains the uptake of N 2 05>. 
Additional simulations focused on the effects of mineral 
dust for the same time period were presented by Martin et 
al. [2002b]. We also perform 1 -month sensitivity simula- 
tions (following a 1 -month spin-up) for March and August 
1997 to compare the relative importance of radiative and 
chemical effects of aerosols. 

3. Photochemical Effects of Aerosols 

[i<>] Figure 3 shows that aerosols reduce UV photolysis 
frequencies near the surface throughout a large part of the 
world. The noontime 0 3 — 0( 1 D) photolysis frequency, 
J(0( 1 D)) (300-320 nrn). decreases by more than 5% over 
most of the Northern Hemisphere, largely due to mineral 



.45 


Figure 3. Sensitivity to aerosols of the 0 3 — O('D) 
photolysis frequency (J(0('D))) in surface air. Values are 
monthly mean model results for March and August 1997. 
See color version of this figure at back of this issue. 
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Figure 4. Fraction of total HO x loss in the lower 
troposphere (~600 m above the surface) contributed by 
uptake of H0 2 by aerosols. Values are monthly mean model 
results for March and August 1997. See color version of this 
figure at back of this issue. 

dust. The effect on J(N0 2 ) (340-400 nni) is similar to that 
on J(0( 1 D)) and about 10% weaker. Removing sulfate from 
the aerosol mixture generally causes a slight increase (<2%) 
in J(0( 1 D)) relative to the standard simulation near the 
surface; 1-D calculations by Liao et al. [1999] similarly 
showed that in the presence of absorbing aerosols, sulfate 
aerosols have a small effect on photolysis frequencies. 
Mineral dust reduces J(0('D)) by 15-30% over and down- 
wind of northern Africa throughout the year [Martin et al., 
2002b], Black carbon causes most of the 15-25% decrease 
in J(0('D)) over northern Europe during August. Aerosols 
from biomass burning and other sources in the Ganges 
Valley reduce J(0( 1 D)) near the surface by a factor of 2, 
largely due to black carbon with smaller contributions from 
organic carbon and mineral dust. Similar features are 
observed over southern Africa during August. The reduc- 
tion in photolysis frequencies may be 5-15% greater for 
internally mixed aerosols than for the external mixture 
assumed here [Liao et al., 1999]. Measurements of 
J(0( 1 D)) over the northern Indian Ocean during March 
showed reductions of 40% due to absorbing aerosols (J. 
Burkert et al., Trace gas and radical diurnal behavior in the 
marine boundary layer during INDOEX 1999, submitted to 
Journal of Geophysical Research , 2002), consistent with the 
results presented here. 

[20] Heterogeneous chemistry also contributes to the 
photochemical effects of aerosols. Figure 4 shows the 
fraction of total HO x loss in the lower troposphere (~600 
m altitude) contributed by aerosol uptake of H0 2 . Here we 
define total HO x loss as the sum of the main loss pathways 
[Frost et al., 1998; Jaegle et al., 2001]: 2(OH + H0 2 ) + 
2(H0 2 + HO.) + 2(H0 2 + CH,0 2 ) + (NO, + OH) + 
2(HN0 4 + OH) + aerosol uptake of H0 2 . Aerosols are 
responsible for more than 70% of HO x loss over biomass 
burning regions. The abundance of fine aerosol over bio- 
mass burning regions [Anderson et al., 1996] and industrial 
regions [ten Brink et al., 1997] leads to high aerosol surface 
area and rapid gas-particle mass transfer (this transfer is not 


diffusion-limited, unlike with supermicron particles). Aero- 
sol uptake of H0 2 is responsible for over 50% of HO x loss 
over Eastern Europe during August, due to high surface 
areas from organic carbon (anthropogenic and biogenic) and 
sulfate, both of which peak in the summer. Black carbon 
constitutes less than 20% of total surface area over Eastern 
Europe, despite playing a dominant role in the reduction of 
photolysis rates. Uptake of H0 2 by aerosols accounts for 
over 20% of HO x loss over much of northern Africa, the 
Atlantic Ocean, and Asia, and over 50% of HO x loss in the 
Arctic during March. 

[21] Figure 5 shows the fraction of total NO x loss to 
HNO3 contributed by aerosol uptake of N0 2 and N0 3 in the 
lower troposphere for March and August. This fraction is 
15-30% in regions with high concentrations of biomass 
burning, mineral dust, and sea-salt aerosols. Over land, 
heterogeneous NO x loss is mainly N0 2 uptake. Over 
oceans, uptake of N0 2 and N0 3 have comparable roles. 
We do not include reaction of NO3 with marine DMS in the 
model, and as a result may exaggerate the importance of 
N0 3 heterogeneous uptake over the oceans. 

[22] We examine the vertical distribution of the photo- 
chemical effects of aerosols over two regions where the 
effects are large: during March over India where carbona- 
ceous aerosols are particularly important and during August 
over the Sahara where mineral dust dominates. The vertical 
profile of the radiative effect is shown in the top panels of 
Figure 6. The effect decreases rapidly with altitude follow- 
ing the distribution of aerosols. The small enhancement in 
J(N0 2 ) in the upper troposphere over India results from 
scattering by aerosols below. The middle panels of Figure 6 
show the vertical profiles of the chemical effect of aerosols. 
Uptake of H0 2 on aerosol surfaces is responsible for about 
80% of HO x loss through the depth of the boundary layer 
over India, and 20-40% of HO x loss throughout most of the 
troposphere over the Sahara. Almost 20% of the HNO3 
production in the boundary layer over India results from 
uptake of N0 2 by aerosols. The bottom panels in Figure 6 
show that over India the radiative and chemical properties 



Figure 5. Fraction of total HNO3 production in the lower 
troposphere (~600 m above the surface) contributed by 
uptake of N0 2 and N0 3 by aerosols. Values are monthly 
mean model results for March and August 1997. See color 
version of this figure at back of this issue. 
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Figure 6. Photochemical effects of aerosols over India and the Sahara. Values are monthly mean model 
results, (top) Effects on J(0('D)) and J(N0 2 ). (middle) Fraction of total HO x loss contributed by aerosol 
uptake of H0 2 and the fraction of total HNO ; production contributed by aerosol uptake of N0 2 and NO,, 
(bottom) Effect of these radiative and heterogeneous processes on OH concentration as determined by 
difference with simulations that did not include either or both of these processes. 


of aerosols have comparable effects on OH near the surface, 
but chemical effects dominate at higher altitudes where the 
overhead aerosol column is smaller. Over the Sahara the 
radiativ e effects of mineral dust are about tw ice as important 
as the chemical effects, reflecting the absorbing nature of 
mineral dust at UV wav elengths [Sokolik et al.. 1993] and 
the large size of the dust particles. 

[ 23 ] Figure 7 shows the combined photochemical effects 
of aerosols on OH. CO. NO x . and O; concentrations in 
March. Here we compare the standard simulation to a 
simulation excluding radiativ e effects of aerosols and the 
heterogeneous reactions ofH0 2 . X0 2 . and NO;, (hydrolysis 
of N 2 0 5 is included in both simulations). The top panel 
shows that OH is reduced by 25-50% over northern Africa 
and South Asia, largely due to mineral dust, organic carbon, 
and black carbon. Over northern India near Calcutta. OH 
decreases by a factor of 4 reflecting comparable contribu- 
tions from radiativ e and chemical effects. Over much of the 
Northern Hemisphere OH concentrations decrease by 5- 


25%: about half of this decrease is from the radiativ e effects 
of mineral dust. The depletion of OH increases the lifetime 
of CO and other trace gases such as NO x and non-methane 
hydrocarbons, increasing their export from regions of emis- 
sions. As a result. CO increases by 5-15 ppbv (5-10%) 
throughout much of the Northern Hemisphere, with local 
increases downwind of biomass burning regions of 20- 
75 ppbv (10—1 5° 0 ). Ship and aircraft measurements during 
the March winter monsoon found elevated CO concentra- 
tions in the marine boundary layer off the coast of India 
[Rhoads et al.. 1997; Leliexeld el al.. 2001]; our analysis 
suggests that depletion of OH by aerosols contributed to the 
observed enhancement. Over the tropical Atlantic down- 
wind of the Sahara. CO increases by 10-20 ppbv (7-13%). 

[24] The photochemical effects of aerosols generally 
result in a slight increase of NO x over continental regions 
due to OH depletion. Over the oceans the effect is more 
often a slight decrease (up to 20%) as NO ; and NO; uptake 
are relatively more important (Figure 5). The bottom panel 
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Figure 7. Changes in OH, CO, NO x , and 0 3 concentrations in the lower troposphere (~600 m above 
the surface) from scattering and absorption by aerosols, and from uptake of HO2, N0 2 , and NO3 by 
aerosols, as determined by difference with a simulation that did not include these effects. Values are 
model monthly means during March. White circles show the locations of observations compared to the 
model in section 4. Also shown is the change in the CO flux due to the above aerosol effects. See color 
version of this figure at back of this issue. 

shows that the effect on O3 is confined largely to regions 
with intense emissions of both aerosols and O3 precursors. 

Ozone is reduced by 10-40% or 15-45 ppbv over and 
downwind of the biomass burning region of northern India 
during March. Ozone also is reduced by 5-10 ppbv (10- 
15%) over northeastern China and the Yellow Sea, a region 
where industrial activity is collocated with high aerosol 
loading from anthropogenic and mineral dust aerosols 
[Zhang et al . , 1994; Zhang and Carmichael , 1999], Phadnis 
and Carmichael [2000] previously used a regional 3-D 


model to calculate a reduction in boundary layer 0 3 of 
about 5% over the region during May due to uptake of 
H0 2 (y = 0.1), H 2 0 2 (y = 0.1), N 2 O s (y = 0.1), and HN0 3 
(y = 0.1) on mineral dust. 

[25] Figure 8 shows the same results but for August. The 
top panel shows that OH is reduced by 25-40% over 
Europe, largely from the radiative and chemical effects of 
sulfate, black carbon, and organic carbon. Similar reduc- 
tions in OH over northern Africa and parts of Asia occur 
primarily from the radiative effects of mineral dust. About 
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Figure 1 . Modeled monthly mean optical depth at 400 nm of different aerosol types for March 1997 
(left column) and August 1997 (right column). Values are calculated from mass concentration fields from 
Ginoux el al [2001] for mineral dust and from Chin et al. [2002] for the other aerosol types. 
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Figure 2. Modeled monthly mean total aerosol surface 
area in the lower troposphere (~600 m above the surface) 
for March and August 1997. Values are calculated from 
mass concentration fields from Ginoux et al. [2001] for 
mineral dust and from Chin et al. [2002] for the other 
aerosol types. 



Figure 3. Sensitivity to aerosols of the 0 3 — *■ 0( 1 D) 
photolysis frequency ( J(0( 1 0))) in surface air. Values are 
monthly mean model results for March and August 1997. 
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Figure 4. Fraction of total HO x loss in the lower 
troposphere (~600 m above the surface) contributed by 
uptake of HCF by aerosols. Values are monthly mean model 
results for March and August 1997. 



Figure 5. Fraction of total HNO, production in the lower 
troposphere (~600 m above the surface) contributed bv 
uptake of NCK and NO;, by aerosols. Values are monthly 
mean model results for March and August 1997. 
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Figure 7. Changes in OH, CO, NO x , and 0 3 concentrations in the lower troposphere (~600 m above 
the surface) from scattering and absorption by aerosols, and from uptake of HO?, NO 2 , and N0 3 by 
aerosols, as determined by difference with a simulation that did not include these effects. Values are 
model monthly means during March. White circles show the locations of observations compared to the 
model in section 4. Also shown is the change in the CO flux due to the above aerosol effects. 
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Figure 8. As in Figure 7. but for August. See color version of this figure at back of this issue. 


half of the reduction of 5-25% over most of the Northern 
Hemisphere is from the radiative effects of mineral dust. 
Concentrations of CO increase bv 5-10 ppbv throughout 
much of the Northern Hemisphere, with local increases of 
more than 10 ppbv over Europe and the Mediterranean. 
Concentrations of NO x generally increase by less than 10%, 
with local increases of up to 30% over Europe, and local 
decreases of up to 20° o over tropical oceans. As for March 
the effect on Ch is localized to photochemically active 
regions where aerosols are collocated with anthropogenic 
emissions of O, precursors. Aerosols reduce 0 3 in the 


boundary layer by 5-9 ppbv over northern Europe during 
August 1997. We find a smaller reduction in July 1997 (3- 
5 ppbv) due to lower aerosol loadings, which in the model 
are due to higher precipitation. The photochemical effects of 
aerosols are smaller over the eastern United States where 
concentrations of absorbing aerosols are lower. The com- 
bined photochemical effects of aerosols decrease 0 3 in the 
boundary layer over the eastern United States by up to 
3 ppbv and over East Asia by up to 5 ppbv. 

[ 26 ] Comparison of our results with previous global 
model studies of sensitivity to aerosol photochemistry' is 
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Table 3. Aerosol Processes Examined in Previous Global Sensitivity Studies 



Aerosol Types 

Included (X) 





Reactions Examined 3 



so% 

BC b 

OC“ 

SS d 

Dust 

7N205 

7H02 

7N02 

7N03 

703 

7HN03 

"yHCHO hv* 

Reference 

X 



X 


0.1 



0.1 




[Dentener and Crutzen , 1993] 





X 

0.1 

0.1 



1-20 x 10~“ 

0.1 


[Dentener et al., 1996] 

X 

X 




0.1 

© 

k> 

1 

p 

c* 



7-300 x 10-“ 


0.01-0.02 

[Tie et at.. 2001] 

X 

X 

X 


X 

0.1 

0.2 

10~ 4 

10' 3 

5 x 10~ 5 

0.1 

X 

[Liao et al 2003] 





X 

0.1 

0.2 

10~ 4 

10~ 3 



X 

[Martin et al 2002a, 2002b] 

X 

X 

X 

X 

X 


0.2 

10- 4 

10' 3 



X 

This work 


“Uptake rates shown only for reactions examined in separate sensitivity simulations. 
b Black carbon. 

“Organic carbon. 
d Sea salt. 

“Scattering and absorption of ultraviolet radiation by aerosols. 


difficult because of differences in the aerosol processes and 
aerosol types examined. Table 3 shows these differences. 
Liao et al. [2003] find that the radiative effects of aerosols 
change 0 3 concentrations by less than 1 ppbv anywhere in a 
global 3-D model with substantially smaller aerosol load- 
ings than used here. Dickerson et al. [1997] calculated that 
the radiative effects of aerosols increase surface 0 3 over 
regions of scattering aerosols such as the eastern United 
States, and decrease 0 3 over regions of absorbing aerosols. 
We find that the combined radiative and chemical effects of 
aerosols lead to small net decreases in 0 3 over the eastern 
United States, and larger decreases in 0 3 over Europe and 
Asia where aerosols are more absorbing. Tie et al. [2001] 
calculated that uptake of H0 2 (y = 0.2) by sulfate aerosols 
decreases zonal mean OH by 5-10% during June north of 
45°N. Over this region sulfate aerosols are the largest 
contributor to the aerosol load; our results are comparable 
to those of Tie et al. [200 1 ] and indicate that consideration 
of radiative effects would about double the impact on OH. 
Dentener et al. [1996] found that uptake of HO, (y = 0.1) 
by mineral dust decreased H0 2 concentrations over north- 
ern Africa and the Gobi Desert by up to 10% during 
February, March, and April, but had little effect on 0 3 
due to low NO x concentrations. Over this region where 
mineral dust dominates the total aerosol load, we calculate 
an effect on H0 2 four times as large using y H 02 of 0.2 and 
accounting for scattering and absorption of UV radiation by 
mineral dust. The effect on 0 3 is still small. 

4. Comparison With Observations 

[27] We compared our model results with a suite of 
aircraft, sonde, and satellite observations to determine if 
inclusion of the photochemical effects of aerosols described 
here improves or degrades the simulation. Again these 
effects include (1) absorption and scattering of ultraviolet 
radiation and (2) reactive uptake of H0 2 , N0 2 , and N0 3 . 
We present here selected comparisons for locations where 
the model exhibits large discrepancies with observations 
and where aerosol effects are large. Large decreases in 0 3 
are found over South Asia in March and northern Europe in 
August due to photochemical effects of aerosols (Figures 7 
and 8). Figure 9 compares monthly mean modeled 0 3 
profiles over South Asia with measurements from commer- 
cial aircraft as part of the Measurement of Ozone and Water 
Vapor by Airbus In-Service Aircraft (MOZAIC) program 
[Marenco et al., 1998] during 1995-1999. The photo- 


chemical effects of aerosols improve the simulation at all 
three sites by decreasing 0 3 concentrations in the lower 
troposphere, although they have minimal effect at higher 
altitudes (above 600 hPa) where the model is still too high. 
The magnitude of the lower tropospheric correction, 
although insufficient to bring model results in agreement 
with observations, is correlated with the magnitude of the 
model bias. The low 0 3 values observed at Madras in 
surface air likely reflect local titration. Lai and Lawrence 
[2001] found a model overestimate of comparable magni- 
tude at the surface during February and March at Ahmeda- 
bad (23°N, 73°E). Our black carbon emissions from fossil 
fuels in India [from Cooke et al., 1999] may be a factor of 
2-3 too low [Dickerson et al., 2002], so that aerosols may 
reduce 0 3 over India throughout the year and may lead to a 
greater 0 3 reduction during March. However, it does not 
appear that aerosols could correct the model overestimate 
above 600 hPa. 

[28] Figure 10 compares simulated and observed (sondes) 
ozone profiles at Hohenpcissenberg in northern Europe 
(43°N, 11°E, 975 m altitude) in August. Topography 
complicates the comparison between GEOS-CHEM and 
ozonesondes near the surface, but the comparison at 800 
hPa shows that aerosols reduce the model bias. There is a 
negative model bias for that site in the free troposphere for 
which aerosols offer no explanation. We calculate that 
aerosols decrease surface 0 3 over the eastern United States 
by 1 -3 ppbv, which could correct our previously reported 
model bias of 3-5 ppbv in the region [Fiore et al., 2002], 

[29] Figure 11 shows that consideration of the photo- 
chemical effects of aerosols increases CO by 2-5 ppbv at 
most tropical sites, by 5-10 ppbv at remote sites in the 
northern midlatitudes, and by about 10 ppbv in the Arctic, 
improving agreement with observations and removing the 
model bias in the extratropics. We find relatively little 
influence on the seasonal variation of CO. 

5. Global Budgets of 0 3 and OH 

[ 30 ] Table 4 shows annual mean global 0 3 production 
and loss rates computed in the model, as well as the lifetime 
of methylchloroform (CH 3 CC1 3 ) against oxidation by tropo- 
spheric OH. The latter is a standard index of the global 
mean tropospheric OH concentration [Spivakovsky et al., 
2000]. The photochemical effects of aerosols beyond N 2 0 5 
hydrolysis reduce global 0 3 production and loss by 6% and 
7%, respectively. The global 0 3 burden is unchanged 
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Figure 9. Comparison of observed and simulated O3 
concentrations over South Asia during March. Lines with 
error bars are aircraft measurements from MOZAIC. 1995- 
1999 [Marenco et al.. 1998], Error bars indicate standard 
deviations. The number of observations is denoted by “n." 
The thick solid line is the monthly mean 0 3 concentration 
from the standard simulation. The dashed line represents 
values from a simulation without the radiative effects or 
uptake of H0 2 , NCK and NO ? by aerosols. 


because the large regional decreases in boundary layer O3 
are compensated by small increases throughout the tropical 
troposphere. The global 0 3 production and loss rates 
calculated here are still at the high end of the current 
generation of global 3-D models, an issue discussed by 
Bey et al. [2001a], Other models, including an earlier 
generation of our model [Wang et al. 1998], indicate global 
production rates in the range 3300-4550 Tg yr~' and 
global chemical loss rates in the range 2500-4100 yr 1 
[ World Meteorological Organization. 1998: Leliexeld and 
Dentener. 2000]. Our higher values relative to Wang et al. 


[1998] reflect more accurate treatment of radiative transfer 
through clouds [Wild et al.. 2000], as discussed by Bey et al. 
[2001a], 

[31] Dentener and Crutzen [1993] previously found that 
N 2 0 5 hydrolysis decreases the global mean tropospheric 
OH concentration (as measured by the lifetime of CH3CCI3) 
by 9%. We find here that the additional photochemical 
effects of aerosols globally decrease OH by an additional 
9%. The radiative effects of mineral dust explain 60% of 
this decrease. The OH decrease is much larger in the 
Northern Hemisphere (13%) than the Southern Hemisphere 
(4%). The resulting global CH3CCI3 lifetime of 5.6 years 
against oxidation by tropospheric OH is consistent with 
estimates from observ ations of 5.7 ± 0.7 years [Spivakovsky 
et al. 2000] and 6.0 (+1.0, —0.7) years [Prim et al. 2001]. 

6. Conclusions 

[3:] We have evaluated the sensitivity of tropospheric 
OH. O3, and O3 precursors to aerosol effects not usually 
included in global models: ( 1 ) aerosol scattering and 
absorption of ultraviolet radiation and (2) reactive uptake 
of HO : , N0 2 . and NO3. We used for this purpose a global 
3-D model of tropospheric chemistry' (GEOS-CHEM) with 
specified global 3-D fields of sulfate, black carbon, organic 
carbon, sea salt, and mineral dust aerosol concentrations 
from a global model (GOCART) driven by the same 
meteorological fields. We did not examine the role of 
N 2 0 5 hydrolysis in aerosols, whose importance is well 
established [Dentener and Crutzen. 1993] and which is 
included in all current-generation tropospheric chemistry 
models. 

[33] We find important aerosol effects from reduction in 
photolysis frequencies and reactive uptake of H0 2 (y = 0.2). 
Perturbations to photolysis frequencies are largely from 
absorbing aerosols: sulfate has little effect. Aerosols 
decrease the O3 — 0< 1 D ) photolysis frequency. J(0( 1 D)), 
at the surface by 5-15% throughout most of the Northern 
Hemisphere, largely due to mineral dust, and by a factor of 
2 in biomass burning regions, largely due to black carbon. 

[34] Uptake of H0 2 by aerosols accounts for more than 
10% of total HO x (= OH + peroxy) loss in the boundary' 
layer of most continental regions. The heterogeneous chem- 
ical effects of aerosols are largest where high concentrations 


Hohenpeissenberg (48N, 11E) 



Ozone (ppbv) 

Figure 10. As in Figure 9. but for August at Hohenpeis- 
senberg (48°N, 11°E, 975 m altitude). Measurements are 
from a 1980-1993 ozonesonde climatology [Logan. 1999], 
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Figure 11 . Comparison of observ ed and simulated annual mean CO concentrations in surface air for 
remote sites (annual mean concentrations less than 140 ppbv) of the Northern Hemisphere. Solid circles 
are observed values for 1988-1996 from Novelli et al. [1998]. Solid line is values from the simulation 
with full aerosol photochemistry. Dashed line is values from a simulation without the radiative effect of 
aerosols and without uptake of H0 2 , N0 2 , and N0 3 by aerosols. 


of fine aerosols enable efficient gas-particle mass transfer. 
Uptake of H0 2 on aerosols accounts for up to 50-70% of 
HO x loss over polluted regions of Eastern Europe (largely 
due to sulfate and organic carbon), and more than 70% over 
biomass burning areas of the tropics (largely due to organic 
carbon). Aerosols decrease boundary layer OH concentra- 
tions by 5-35% in most of the Northern Hemisphere, by a 
factor of 2 over northern Europe during August, and up to a 
factor of 4 over India during March. 

[ 35 ] Reactive uptake of N0 2 (y = 10“ 4 ) and N0 3 (y = 
10’) have relatively small effects. These processes account 
for more than 10% of HN0 3 production only over the 
tropical North Atlantic, the Sahara, and the southern oceans. 
Over continental regions, the net effect of the aerosol 
processes considered here is to slightly increase NO x due 
to OH depletion. 

Table 4. Global Tropospheric 0 3 Budget and CH 3 CCI 3 Lifetime 2 

Full Aerosol N 2 0 5 Hydrolysis 
Photochemistry Only 11 


Chemical production, Tg O, yr 

4924 

5263 

Chemical loss, Tg O, yr 

4377 

4693 

Burden, Tg O 3 

322 

323 

Global CH 3 CCI 3 lifetime, yr 

5.6 

5.2 

NH CH,CC1, lifetime, yr 

5.4 

4.8 

SH CH 3 CCI 3 lifetime, yr 

5.8 

5.6 


“The production, loss, and burden of O, are actually for the extended odd 
oxygen family defined as 0 3 + N0 2 + 2 * NO, + peroxyacylnitrates + 
HNO, + 3 * N 2 O s . They are calculated for the column extending up to the 
local model tropopause. The CH,CC1, lifetime is calculated as the ratio of 
the total burden of atmospheric CH,CC1 3 to the tropospheric loss rate 
against oxidation by OH [Spivakovsky et al, 2000]. Values are annual 
means for September 1996 to August 1997 from the GEOS-CHEM model. 

b This is the same as full aerosol photochemistry simulation but without 
radiative effects or uptake of NO : . NO,, and HO. by aerosols. 


[36] We examined the implications of the above aerosol 
effects on the global model budgets of OH, 0 3 , and CO. 
Annual mean OH concentrations decrease by 9% globally 
(comparable to the effect of N 2 0 5 hydrolysis) and by 13% 
in the Northern Hemisphere. About 60% of the global 
decrease is due to the radiative effects of mineral dust. 
The resulting CH 3 CC1 3 lifetime against oxidation by tropo- 
spheric OH (5.6 years) is in better agreement with estimates 
constrained by observations. Annual mean global 0 3 chem- 
ical production decreases by 6%, but the tropospheric 0 3 
burden remains unchanged. Concentrations of CO increase 
by 5- 15 ppbv in most the Northern Hemisphere, improving 
agreement with observations at remote sites. Boundary layer 
0 3 decreases by 15-45 ppbv over South Asia during the 
biomass burning season in March, improving the compar- 
ison to aircraft observations (MOZAIC) although the model 
is still too high. Summertime boundary layer 0 3 is reduced 
by 5-9 ppbv over northern Europe, and by 1-3 ppbv over 
the United States, again improving the model simulation 
with respect to observations over both regions. We conclude 
that surface 0 3 concentrations over Europe and other 
industrial regions will increase if emissions of aerosols are 
reduced without corresponding reductions in 0 3 precursors. 
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